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Abstract 

We argue that the critical behaviour near the point of "gradient catastrophe" of 
the solution to the Cauchy problem for the focusing nonlinear Schrodinger equation 
it ^ t + y^xa; + l^l 2 ^ = 0, e « 1, with analytic initial data of the form ^(x, 0; e) = 
A(x) e? 5 ^' is approximately described by a particular solution to the Painleve-I equa- 
tion. 

1 Introduction 

The focusing nonlinear Schrodinger (NLS) equation for the complex valued function ^ = 

t q t + ^ xx + \q\^ = (1.1) 

has numerous physical applications in the description of nonlinear waves (see, e.g., the books 
ll47l l35l l36l ). It can be considered as an infinite dimensional analogue of a completely 
integrable Hamiltonian system [|48l . where the Hamiltonian and the Poisson bracket is given 
by 

% + {V(x),H} = 

{V(x),**(y)} = i5(x-y) (1.2) 

n = \J (|*,| 2 -|*l 4 ) dx 



i 



(here ty* stands for the complex conjugate function). Properties of various classes of so- 
lutions to this equation have been extensively studied both analytically and numerically 
JIllSCZllSllI3EIlEIl|2gl^lll2llSlB3E|. One of the striking features that distinguishes 
this equation from, say, the defocusing case 



-\fj 



\m 2 m = o 



is the phenomenon of modulation instability 03 [EH HH. Namely, slow modulations of the 
plane wave solutions 



develop fast oscillations in finite time. 

The appropriate mathematical framework for studying these phenomena is the theory of 
the initial value problem 

tf(x,0;e) = A(x)e" s{x) (1.3) 

for the e-dependent NLS 

ie% + -* lgx + \*\ 2 * = 0. (1.4) 

Here e > is a small parameter, A(x) and S(x) are real-valued smooth functions. Introduc- 
ing the slow variables 



u 



2i V * 

the equation can be recast into the following system: 

u t + (uv) x = 



X 



(1.5) 



(1.6) 



V t + v V x - u x + 



4 \ 2m 2 



U 



0. 



The initial data for the system ( |1.6[ ) coming from ( |1.3[ ) do not depend on e: 

u(x, 0) = A 2 (x), v(x, 0) = S'(x). 



(1.7) 



The simplest explanation of the modulation instability then comes from considering the so- 
called dispersionless limit e — > 0. In this limit one obtains the following first order quasilin- 
ear system 

u t +v u x + u v x — 



(1.8) 







This is a system of elliptic type because of the condition u > 0. Indeed, the eigenvalues of 
the coefficient matrix 



v u 
-1 V 
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are complex conjugate, A = v ± iy/u. So, the Cauchy problem for the system ( |1.8| ) is ill- 
posed in the Hadamard sense (cf. |33l|7l). Even for analytic initial data the life span of a 
typical solution is finite, t < t . The x- and t-derivatives explode at some point x = xq when 
the time approaches to- This phenomenon is similar to the gradient catastrophe of solutions 
to nonlinear hyperbolic PDEs [2 J . 



For the full system (1.6) the Cauchy problem is well-posed for a suitable class of e- 



independent initial data (see details in |[T7ll46lO . However, the well-posedness is not uniform 



in e. In practical terms that means that the solution to ( |1.6[ ) behaves in a very irregular way 
in some regions of the (x,t)-plane when e — > 0. Such an irregular behaviour begins near 
the points (x = xq, t = to) of the "gradient catastrophe" of the solution to the dispersionless 



limit ( |1.8[ ). The solutions to ( | 1 . 8[ ) and ( |1.6[ ) are essentially indistinguishable for t < to', the 
situation changes dramatically near x when approaching the critical point. Namely, when 
approaching t = t the peak near a local maximum^ of u becomes more and more narrow 
due to self-focusing; the solution develops a zone of rapid oscillations for t > t . They 
have been studied both analytically and numerically in [f8l fT6l I2T1 [23l l25l [26l l34l l43l l44i 
However, no results are available so far about the behaviour of the solutions to the focusing 
NLS at the critical point (xo, t ). 

The main subject of this work is the study of the behaviour of solutions to the Cauchy 
problem ( |1.6[ ), ( |1.7[ ) near the point of gradient catastrophe of the dispersionless system ( |1.8[ ). 



In order to deal with the Cauchy problem for < | 1 . 8[ ) we will assume analytichyj^of the initial 



data u(x, 0), v (x, 0). Then the Cauchy problem for ( |1.8[ ) can be solved for t < t via a 
suitable version of the hodograph transform (see Section [2] below). An important feature of 
the gradient catastrophe for this system is that it happens at an isolated point of the (x,t)- 
plane, unlikely the case of KdV or defocusing NLS where the singularity of the hodograph 



solution takes place on a curve. We identify this singularity for a generic solution to ( |1.8[ ) 
as the elliptic umbilic singularity (see Section [4] below) in the terminology of R.Thom ||42~1|. 
This codimension 2 singularity is one of the real forms labeled by the root system of the D± 
type in the terminology of V.Arnold et al. J3J. 

Our main goal is to find a replacement for the elliptic umbilic singularity when the disper- 
sive terms are added, i.e., we want to describe the leading term of the asymptotic behaviour 



for e — > of the solution to ( |1.6[ ) near the critical point (x , t ) of a generic solution to ( |1.8[ ). 

Thus, our study can be considered as a continuation of the programme initiated in IfTBl 
to study critical behaviour of Hamiltonian perturbations of nonlinear hyperbolic PDEs; the 



fundamental difference is that the non perturbed system < | 1 . 8| ) is not hyperbolic! However, 
many ideas and methods of lfT3l (see also [12]) play an important role in our considerations. 

The most important of these is the idea of universality of the critical behaviour. The 



'Regarding initial data with local minima we did not observe cusps related to minima in numerical simula- 
tions. We believe they do not exist because of the focusing effect in the NLS that pushes maxima to cusps but 
seems to smoothen minima. 

2 We believe that the main conclusions of this paper must hold true also for non analytic initial data; the 
numerical experiments of [ 8 1 do not show much difference in the properties of solution between analytic and 
non analytic cases. However, the precise formulation of our Main Conjecture has to be refined in the non 
analytic case. 



3 



general formulation of the universality suggested in |fT3l for the case of Hamiltonian pertur- 
bations of the scalar nonlinear transport equation says that the leading term of the multiscale 
asymptotics of the generic solution near the critical point does not depend on the choice of 
the solution, modulo Galilean transformations and rescalings. This leading term was identi- 
fied in lfT3l via a particular solution to the fourth order analogue of the Painleve-I equation 
(the so-called Pf equation). The existence of the needed solution to Pf has been rigorously 
established in [9|. Moreover, it was argued in [13] that this behaviour is essentially indepen- 
dent on the choice of the Hamiltonian perturbation. Some of these universality conjectures 
have been partially confirmed by numerical analysis carried out in 11211 . 

The main message of this paper is the formulation of the Universality Conjecture for 
the critical behaviour of the solutions to the focusing NLS. Our considerations suggest the 



description of the leading term in the asymptotic expansion of the solution to ( |1.6[ ) near the 
critical point via a particular solution to the classical Painleve-I equation (P-I) 

n C( = 6 ft 2 -c 

The so-called tritronquee solution to P-I was discovered by PBoutroux as the unique 
solution having no poles in the sector | arg£| < 4.7r/5 for sufficiently large \C\. Remarkably, 
the very same solution^] arises in the critical behaviour of solutions to focusing NLS! 

The paper is organized as follows. In Section [2] we develop a version of the hodograph 
transform for integrating the dispersionless system ( |1.8[ ) before the catastrophe t < to. We 



also establish the shape of the singularity of the solution near the critical point; the latter is 
identified in Section [4] with the elliptic umbilic catastrophe of Thorn. In Section|3]we develop 
a method of constructing formal perturbative solutions to the full system ( |1.6[ ) before the 
critical time. In Section|5]we collect the necessary information about the tritronquee solution 
of P-I and formulate the Main Conjecture of this paper. Such a formulation relies on a much 
stronger property of the tritronquee solution: namely, we need this solution to be pole-free 
in the whole sector | arg£| < Aix/h. Numerical evidence for the absence of poles in this 
sector is given in Section [6} In Section [7] we analyze numerically the agreement between the 
critical behaviour of solutions to focusing NLS and its conjectural description in terms of the 
tritronquee solution restricted on certain lines in the complex £-plane. In the final Section [8] 
we give some additional remark and outline the programme of future research. 
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2 Dispersionless NLS, its solutions and critical behaviour 



The equations ( |1.6[ ) are a Hamiltonian system 

ut + {u(x),H} = 

v t + {v{x),H} = 
with respect to the Poisson bracket originated in ( |1.2[ ) 

{u(x),v(y)} = 8'(x-y), 
other brackets vanish, with the Hamiltonian 



(2.1) 



H 



1 



(uv 2 -u 2 )+ c 2 



8u 



dx. 



(2.2) 



Let us first describe the general analytic solution to the dispersionless system ( 1.8 ). 



Lemma 2.1 Let uj(x), vj(x) be two real valued analytic functions of the real variable x 
satisfying 

(Ul,x) 2 + {Vl,xf ^ 0. 

Then the solution u = u(x, t), v = v(x, t) to the Cauchy problem 

u(x, 0) = uj(x), v (x, 0) = vi(x) (2.3) 
for the system { 1.8\ for sufficiently small t can be determined from the following system 

x = vt + f u 

= ut + f v 

where f = f(u, v) is an analytic solution to the following linear PDE: 

fw + uf uu = 0. (2.5) 



(2.4) 



Conversely, given any solution to ( 2.5) sa tisfying u 2 f 2 u + f 2 v ^ at some point (u = u,v = 
v) such that f v (u, v) = 0, the system (2.4) determines a solution to (1.8) defined locally near 



the point x = x := f u (u, v)for sufficiently small t. 
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Remark 2.2 The solutions to the linear PDE ( 2.5 ) correspond to the first integrals ofdis- 
persionless NLS: 

F = J f(u,v)dx, j f F = 0. (2.6) 
Taking them as the Hamiltonians 

u s + {u(x),F} = u t + (f v ) x = 

v s + {v(x),F} = v t + (f u ) x = 
yields infinitesimal symmetries of the dispersionless NLS: 

(u t ) s = (u s ) t , (v t ) s = (v s ) t . 



(2.7) 



(2.8) 



One of the first integrals will be extensively used in this paper. It corresponds to the 
Hamiltonian density 

(2.9) 



9 = ~v 2 + u(\ogu - 1). 



The associated Hamiltonian flow reads 



u + v — 



(2.10) 



Eliminating the dependent variable v one arrives at the elliptic version of the long wave limit 
of Toda lattice: 

u ss + (\ogu) xx = 0. 

Due to commutativity (2.8) the systems (1.8) and ( 2.10[ ) admit a simultaneous solution u = 
u(x, t, s), v = v(x, t, s). Any such solution can be locally determined from a system similar 
to d24]) 

X = V t + f u 

' (2.11) 

s = u t + f v 



where / = f(u, v), as above, solves the linear PDE (2.5 ). 



The system (2.1 1 ) determines a solution u = u(x, t, s), v = v(x, t, s) provided applica- 
bility of the implicit function theorem. The conditions of the latter fail to hold at the critical 
point (xq, so,t ,uo,vo) such that 



^0 = ^0*0 + fu(u ,v ) 
so = u t + f v (uo,v ) 

fuu(u ,v ) = f vv {u ,v ) = 0, f uv (u ,v ) = -t 



(2.12) 
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In sequel we adopt the following system of notations: the values of the function / and of its 
derivatives at the critical point will be denoted by f° etc. E.g., the last line of the conditions 



(2.12) will read 



f = f = n f° = —t 

Definition 1. We say that the critical point is generic if at this point: 

f° ± 0. 

J UUV I " 

Let us the introduce real parameters r, ip determined by the third derivatives of the function 
/ evaluated at the critical point, 



(cos ip - i sin ip) = f° uv + i^/uof, 



uuw 



(2.13) 



Due to the genericity assumption 



7T 



In order to describe the local behaviour of a solution to the dispersionless NLS/Toda 
equations we define a function R(X;S,ip) of real variables X, S depending on the real 
parameter ip satisfying 

(5 + cos^) 2 + (X + sinV)V ( 2 - 14 ) 



by the following formula 

R(X, S, ip) = sign [cos ip] 

xJl + X simp + S costp + ^1 + 2(X sint/j + S cosi(j) + X 2 + S 2 



(2.15) 



Put 



Qo(X,S^) 



V2 



V2 



R(X, S,i[)) cos ip 



(X cos ip — S sin ip) sin ip 



(X cos ip — S sin ip) cos ip 

R(x,S,ip) 



R(X,S,iP) 



+ R(X, S,ip) sinip 



cos ip 



(2.16) 



sinip. 



Observe that Po(X, S, ip) and Qo(X, S, ip) are smooth functions of the real variable X pro- 
vided validity of the inequality ( |2.14[ ). 



Lemma 2.3 Given an analytic solution u(x, s, t), v(x, s, t) to the dispersionless NLS/Toda 
equations with a generic critical point (x , Sq, t , Uq, Vq), and arbitrary real numbers X, S 
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satisfying $2.14\, T < 0, then there exist the following limits 



lim \- 1 ' 2 

A^+O 



a u + X 1/2 v T + ^=rXT 2 , s + X 1/2 u T + ^-rST 2 , A 1/2 T) - u 



rTP (X,S,ij) 



lim \- 1 ' 2 

A^+0 



(2.17) 

r(,: + X 1/2 v T + ^rXT 2 , s + X 1/2 u T + ^-rST 2 , X 1/2 T) - v 



TQ Q (X, 



where the parameters r, ip are defined by (2.13). 



Proof From the linear PDE (2.5 ) it follows that 



fu 



'llfuuu fuut fvvv ufu 



Using these formulae we expand the implicit function equations ( |2.1 1\ near the critical point 
in the form 



x 



~ v i=vi+ \ [f uuu (u 2 - u v 2 ) + 2f uuv uv] +0((\u\ 2 + \v\y/ 2 ) 



s - u i=ui+ \ [f uuv (u 2 - u v 2 ) - 2u f uuu uv] +0{(\u\ 2 + \v\ 2 f 2 ) 



(2.18) 



where we introduce the shifted variables 



The rescaling 



X = X — Xo, S = S — So, t = t — to 
U = U — Uq, V = V — Vq. 



X — v t 

S — Uot 

i 
u 



X(x — v i) 
X(s — u t) 
X^H 
X 1/2 u 



(2.19) 



v h- > X l l 2 v 

in the limit A — ► yields the quadratic equation 

1 2 _ . 

z = tw + -aw , t 7& 



(2.20) 
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where the complex independent and dependent variables z and w read 

Z = S + 1\/UqX — (llQ + ^^/Uo^>o)t, W = U + iy/U()V 

and the complex constant a is defined by 

a — fuuv + iV^ofuum 

therefore 



(2.21) 



(2.22) 



1 



r e 



The substitution 



rt 2 



rt 2 



reduces the quadratic equation to 

(w + ^e^) 2 =r¥e^ [l + e^* (S + i X 
For i < we choose the following root 



w = rte 



y/l + e-^iS + iX) - 1 



(2.23) 



where the branch of the square root is obtained by the analytic continuation of the one taking 
positive values on the positive real axis. Equivalently, 



w = tre** 
where 



sign (cos if)) — ( a/A + 1 + S cos if) + X sin if) + i 
v2 V 



X cos if) — S sin if) 



V A + 1 + S cos if; + X sin if) 



A = a/1 + 2(S'cos^ + XsinV) + X 2 + S 2 . 



This gives the formulae (2.15 ) 



□ 



The result of the lemma describes the local structure of generic solutions to the dis- 
persionless NLS/Toda equations near the critical point. It can also be represented in the 
following form 



u(x,s,t) ~u + rTP (X, S,ip) 



v(x, s,t)~v + —rT Q (X, S, if)) 



(2.24) 



where 



X = 2v ^^, S = 2 £^~ T = t. 



rt 2 



rt 2 ' 



(2.25) 
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We want to emphasize that the approximation ( |2.24[ ) works only near the critical point 



Indeed, for large x — > ±00 the function u(x,s,t) and v(x,s,t) have the following behaviour 



u = — \ r \x it 



'r \x 



\I A a/1 =F shv0 + Mo — r £ cos ^ + O ( . ) (2.26) 

VvMy 

ign (cos V>) \/i ± sin ip +v T — i sin ip + O \ \ ] . (2.27) 

V u o \vW\J 



So, for sufficiently large \x\ the function u(x, s, t) defined by ( |2.24[ ) becomes negative. 
The function u has a maximum at the point X = S tanip, so locally 



u < u + r Tcos ip — \fr I cos-^l \ 1 2 + r T 2 < it - (2.28) 

cos ^ 

At the critical point (x , s , tcb M (b w o) the function u develops a cusp. Let us consider only 
the particular case S = in order to avoid complicated expressions. In this case the local 
behaviour of the function u near the critical point is given by 



u — yjr \x\ y/1 — sirr0, X > 
lim u = < (2.29) 

t-4-Q 



u o ~ \A \ x \ VI + sin?/', x < 

(here £ = y/u Q {x — ^o^))- Thus the parameters r, t/> describe the shape of the cusp at the 
critical point. 



3 First integrals and solutions of the NLS/Toda equations 



Let us first show that any first integral (2.6) of the dispersionless equations can be uniquely 
extended to a first integral of the full equations. 



Lemma 3.1 Given a solution f = f(u, v) to the linear PDE (2.5), there exists a unique, up 
to a total derivative, formal power series in e 

h f = f + J2 ^hf («, v- u x , v x ,..., u^ k \v^) 
fc>i 

such that the integral 

Hj = J hf dx 
commutes with the Hamiltonian ofNLS equation: 

{H,H f } = 



10 



at every order in e. Explicitly, 



hf = f 



,2 r 



12 



(3.1) 



1 

120 



fuuuu + ^ fuuu ) U xx + 2f mUxxVxx U f uuuu V 



1 



fuUUuU XX V x fuUUvV XX U x . o (30 fuUU 9M fuuUU 12M f§ u ~\~ All f(ju) 

80 48m SAbbu 6 



u ; 



1 



432m 2 
1 



f ^fuuuv ^uf uuuuv -\- 2tl f§uv) ^ x ^x 00 o (9 fuuuu ^u/jjtt -\- 2tl fou) ^£^2 

ZOOM 



^~oicn i^fuuuuv ~t~ lO^/kwu 

) u z i£ - — — (18/ 5u + 5m/ 6u ) v x \ + 0(e b ) 



2160 



4320 



Here we use short notations 



fhu ■ 



du 5 ' 



6u 



fl 



5u v 



d'f 

du 5 dv 



Example 1. Taking / = \(uv 2 — u 2 ) one obtains the Hamiltonian of the NLS equation 

hf = \{uv 2 -u 2 )+ e -ul. 

In this case the infinite series truncates. It is easy to see that the series in e truncates if 
and only if f(u, v) is a polynomial in u. Polynomial in u solutions to the linear PDE (2.5 ) 
correspond to the standard first integrals of the NLS hierarchy. 



Example 2. Taking g = —\v 2 + w(logw — 1) (cf. (2.9)) one obtains the Hamiltonian of 
Toda equation 

1 e 2 
h g = --v 2 + u(\ogu - 1) - ^ (u 2 x + 2m v 2 x ) 



-A- 



2 2 2 

U XX _|_ V xx _|_ U xx v x 



+ 



\2A0u 3 60m 2 40m 3 144u 5 24m 4 360m 3 
written in terms of the function <\> = log u in the form 

e 2 <J)xx + e 0(s+e) - 2e^ s) + e^ s ~ e) = 0. 



+ 0(e fc 



(3.2) 



Lemma 3.2 Any solution to the NLS/Toda equations in the class of formal power series in e 
can be obtained from the equations 

6H f 



x = v t + 



s = ut 



5u(x) 
Sv(x) 



(3.3) 
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where f = f(u,v;e) is an arbitrary admissible solution to the linear PDE (2.5) in the class 
of formal power series in e, 

Hf = I hf dx. 



Now, we can apply to the system (3.3 ) the rescaling (2.19) accompanied by the transfor- 
mation 

e i-> A 5/4 e. (3.4) 
At the limit A — > we arrive at the following system of equations 

s - u t = u t + f° u 



\(U 2 - UqV 2 ) + e —U xx 



uof, 





uuu 



uv+ —v xx 





(3.5) 



x — V t = V t + f 





uuu 



yU 2 ~ U V) + —U xx 

I o 



f° 

J uu 



uv+ —v xx 
o 



Using the complex variables z, w defined in (2.21 ) we can rewrite the system in the following 
form: 



z re 



i/U> 



wtre^ + ^w 2 + —w xx . 



(3.6) 



1 

2 U ' 6 

The last observation is that the Toda equations generated by the Hamiltonian H g — J h g dx 
(see Example 2 above) after the scaling limit (2.19 ), ( 3.4| ) yield the Cauchy - Riemann equa- 
tions for the function w = w(z), 

dw/dz = 0. 

Therefore the system ( |3.5[ ) can be recast into the form equivalent to the Painleve-I (P-I) 



equation (see (5.1 ) below) 



z re 



iip 



wtre 1 ^ + -w 
2 

A = e 4 / 5 



Choosing 



l 2 e 



U W ZZ . 



(3.7) 



we eliminate e from the equation. 

In the Section [5] below we will write explicitly the reduction of ( |3.7[ ) to the Painleve-I 
equation and give a conjectural characterization of the particular solution of the latter. 



4 Critical behaviour and elliptic umbilic catastrophe 



Separating again the real and complex parts of (3.6) one obtains a system of ODEs 



e 2 1 

—Uxx + -(U 2 - V 2 ) +rt (Ucos^ - V sinip) -r (S cosip - X sinip) = 
6 2 



(4.1) 



6 



Vxx + UV + r t (U sin if) + V cos if)) — r (S sin if) + X cos if)) = 



12 



that can be identified with the Euler - Lagrange equations 



SS = 0, S = J C(U, V, U x , V x ) dX 

with the Lagrangian 

u x) + \ {U 3 -3UV 2 ) + \rt \(U 2 - V 2 )cosif) -2UVsmif\ 
6 2 

(4.2) 

S cos ip)U + r (S sin if) + X cos if)) V. 

In the "dispersionless limit" e — >• the Euler - Lagrange equations reduce to the search of 
stationary points of a function (let us also set t = 0) 

F =- (U 3 -3UV 2 ) +a + U + a^V (4.3) 
6 

where we redenote 

a + = r (X sin t/> — S cos ■?/>) , a_ = r (S 1 sin if) + X cos ■?/>) . 

At a + = a_ = the function F has an isolated singularity at the origin U = V = of the 
type D 4 also called elliptic umbilic singularity, according to R.Thom [|42l . This singular- 
ity appears in various physical problems; we mention here the caustics in the collisionless 
dark matter 11401 to give just an example. The parameters a + and a define two particular 
directions on the base of the miniversal unfolding of the elliptic umbilic; the full unfolding 
depending on 4 parameters reads 

F = - (U 3 - 3U V 2 ) + -b{U 2 + V 2 ) + a+U + a^V + c. (4.4) 
6 2 

It would be interesting to study the properties of the modified Euler - Lagrange equations for 
the Lagrangian 

c = c + h(u 2 + v 2 ). 

This deformation does not seem to arrive from considering solutions to the NLS hierarchy. 

5 The tritronquee solution to the Painleve-I equation and 
the Main Conjecture 

In this section we will select a particular solution to the Painleve-I (P-I) equation 

n cc = m 2 -( (5.i) 




+r (X sinif) — 
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Recall [1221 that an arbitrary solution to this equation is a meromorphic function on the com- 
plex £-plane. According to P. Boutroux flU the poles of the solutions accumulate along the 
rays 

ar gC= , n = 0, ±1, ±2. (5.2) 



Boutroux proved that, for each ray there is a one-parameter family of particular solutions 
called integrates tronquees whose lines of poles truncate for large £. He proved that the 
integrate tronquee has no poles for large |£| within two consecutive sectors of the angle 
2tt/5 around the ray, and, moreover it has the asymptotic behaviour of the form 



Q = - 



1/2 



l + O 



(5.3) 



for a suitable choice of the branch of the square root (see below) and a sufficiently small 

e > 0. 



Furthemore, if a solution truncates along any two of the rays ( |5.2[ ) then it truncates along 
three of them. These particular solutions to P-I are called tritronquees. They have no poles 



for large |£| in four consecutive sectors; their asymptotic s for large £ is given by ( |5.3[ ). It 
suffices to know the tritronquee solution f2o(C) f° r tne sector 



argCI < 



47T 



(5.4) 



In this case the branch of the square root in ( |5.3[ ) is obtained by the analytic continuation 
of the principal branch taking positive values on the positive half axis £ > 0. Other four 
tritronquees solutions are obtained by applying the symmetry 



0.(0 



4 7T i n 

e s 



2 n I It 

e s 



n 



±1, ±2. 



(5.5) 



The properties of the tritronquees solutions in the finite part of the complex plane were 
studied in the important paper of N.Joshi and A.Kitaev [|24l . 



A. Kapaev [27J obtained a complete characterization of the tritronquees solutions in 
terms of the Riemann - Hilbert problem associated with P-I. We will briefly sketch here 
the main steps of his construction. 



The equation ([5TT]) can be represented as the compatibility condition of the following 
system of linear differential equations for a two-component vector valued function $ = 

*(A,C) 



n ( 2A 2 + 2fiA-C + 2fi 2 



2(x-n) -n c 



(5.6) 



X + 2Q 



(5.7) 
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The canonical matrix solutions $fc(A, C) to the system ( |5.6[ ) - ( |5.7[ ) are uniquely determined 
by their asymptotic behaviour 



S*(A,C) 



1 



A l/4 



V A- 1/4 



A V4 
_ A -l/4 



1 



1 



in the sectors 
Here 



1 

2A 



A G C 



# 2 O 

o h 2 



2tt 



O (A~ 3/2 ) 











(5.8) 



oo, A G Sfc 



^(A,C) = -A 5 / 2 -CA 1 



/2 



0"3 



< are A < 



1 
-1 



2tt 



k + 



keZ. 



(5.9) 



-O 2 
2 c 



2 3 + CO, (5.10) 



the branch cut on the complex A-plane for the fractional powers of A is chosen along the 
negative real half-line. 

The Stokes matrices Sk are defined by 



$ fc+1 (A, C) = $ fc (A, C) S k , A G S fe n E fc+ i. 
They have the triangular form 



5- 



2fc-l 



1 «2fc-l 

1 



S- 



2k 



1 



and satisfy the constraints 

Sk+5 = 0\ Sk 

where 



0"i 



A; G Z; S1S2S3S4S5 = io\ 



1 

1 



(5.11) 



(5.12) 



(5.13) 



Due to ( |5.13| ) two of the Stokes multipliers Sk determine all others; they depend neither on A 
nor on £ provided 0(C) satisfies ( |5.1[ ). 



In order to obtain a parametrization of solutions to the P-I equation ( |5.1[ ) by Stokes mul- 
tipliers of the linear differential equation ( |5.6| ) one has to reformulate the above definitions 
as a certain Riemann - Hilbert problem. The solution of the Riemann - Hilbert problem de- 
pends on C through the asymptotics ( |5.8| ). If the Riemann - Hilbert problem has a unique 
solution for the given Co G C then the canonical matrices $fc(A, C) depend analytically on 
C for sufficiently small |C — Co 1 5 me coefficient O = 0(C) will then satisfy ( |5.1[ ). The poles 
of the meromorphic function 0(C) correspond to the forbidden values of the parameter C for 
which the Riemann - Hilbert problem admits no solution. 
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We will now consider a particular solution to the P-I equation specified by the following 
Riemann - Hilbert problem. Denote four oriented rays 70, 7±i, p in the complex A-plane 
defined by 



7 fc = {A G C I arg A 



2nk 



}, k = 0,±l 



(5.14) 



p = {A G C I arg A = 7r} 



directed towards infinity. The rays divide the complex plane in four sectors. We are looking 
for a piecewise analytic function $(A, () on 

A G C \ ( 7 _! U 70 U 7l U p) 

depending on the parameter C, continuous up to the boundary with the asymptotic behaviour 
at I A I — > 00 of the form ( |5.8[ ) satisfying the following jump conditions on the rays 

$ + (A,C) = <t>-(A,C)<Sfc, \e lk 



(5.15) 



$ + (A,C) = MA,C)S'p, Xe P- 



Here the plus/minus subscripts refer to the boundary values of $ respectively on the left/right 
sides of the corresponding oriented ray, the jump matrices are given by 



1 

% 1 



1 i 
1 



S n 



-i 
-i 



(5.16) 



The following result is due to A.Kapae\|^]. 



Theorem 5.1 The solution to the above Riemann - Hilbert problem exists and it is unique 
for 

47T 



are;A| < 



|A| > R 



(5.17) 



for a sufficiently large positive number R. The associated function 
dH(() 



fio(C) 



H(0 



(5.18) 



11 



is analytic in the domain ( |5.17[ ), it satisfies P-I and enjoys the asymptotic behaviour 

Air 



no(C) 



ICI 



00, 



are A < 



(5.19) 



Moreover, any solution of P-I having no poles in the sector ( |5.17[ ) for some large R > 
coincides with O (C)- 



4 Our solution flo(C) coincides with yz(x) = y-2(x), x = of l27l (see eq. (2.73) of [27 1 ; Kapaev uses 
a different normalization y" — 6y 2 + x of the P-I equation). 
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Joshi and Kitaev proved that the tritronquee solution has no poles on the positive real 
axis. They found a numerical estimate for the position of the first pole Co of the tritronquee 
solution Oo(C) on the negative real axis: 

Co ^ -2.3841687 

(cf. also IfTOl ). Besides this estimate very little is known about the location of poles of fio(0- 
Our numerical experiments (see below) suggest the following 

Main Conjecture. Part 1. The tritronquee solution VLq{Q has no poles in the sector 

Air 

argA < — . (5.20) 
5 



We are now ready to describe the conjectural universal structure behind the critical be- 
haviour of generic solutions to the focusing NLS. For simplicity of the formulation let us 
assume cosip > 0. 

Main Conjecture. Part 2. Any generic solution to the NLS/Toda equations near the 
critical point behaves as follows 

u(x, s,t; e) + iy/u^v(x, s,t; e) ~ u + iy/uoV Q - tre** + 2 e 2/5 (3r v /^) 2/5 e^fi (C) + O (e 

(5 



3r 



1/5 



u t + i^/u^(x — v t) + ^re^t 2 

^75 



where ^o(C) zs the tritronquee solution to the Painleve -I equation (5.1 ). 

The above considerations can actually be applied replacing the NLS time flow by any 
other flow of the NLS/Toda hierarchy. The local description of the critical behaviour remains 
unchanged. 

Remark 5.2 Note that the angle of the line Q(x) in ( 5.21\ for i fixed is equal to ijj/5 + n/2, 
ijj G [— 7T, 7r]. Thus the maximal value ofarg( is equal to 77r/10 < Ati/5. The lines in (5.21 ) 
consequently do not get close to the critical lines of the tritronquee solution. 



6 Numerical analysis of the tritronquee solution of P-I 

In this section we will numerically construct the tritronquee solution VLq, i.e. the tritronquee 



solution with asymptotic behavior ( |5.3[ ). We will drop the index in the following. The 
solution will be first constructed on a straight line in the complex plane. In a second step 
we will then explore global properties of these solutions within the limitations imposed by a 
numerical approacrj^J 

5 Cf. [ 15 1 where a similar technique was applied to solve numerically the Painleve-II equation in the complex 
domain. 
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Let the straight line in the complex plane be given by £ = ay + b with a, b G C constant 
(we choose a to have a non-negative imaginary part) and i/6l. The asymptotic conditions 
are 

for y — > ±00. The root is defined to have its cut along the negative real axis and to as- 
sume positive values on the positive real axis. This choice of the root implies the following 
symmetry for the solution: 

Q(C*)=0*(C). (6.2) 

Thus f2 is real on the real axis, see ll24"ll . 

Numerically it is not convenient to impose boundary conditions at infinity. We thus 
assume that the wanted solution can be expanded in a Laurent series in around infinity. 
Such an asymptotic expansion is possible for the considered tritronquee solution in the sector 
I arg (\ < Air/ 5. The formal series can be written there (see [24]) in the form 



nf = -\iY.T£x, (6-3) 




where a = 1, and where the remaining coefficients follow from the recurrence relation for 

k > 

25k 2 -11* 

a k+l — 7= O-k — 77 > ^mO-k+l-m- (6.4) 

8V6 * — \ 

m=l 

This formal series is divergent, the coefficients behave asymptotically as ((k — l)!) 2 , see 
[1241 for a detailed discussion. 

It is known that divergent series can be used to get numerically acceptable approxima- 
tions to the values of the sum by properly truncating the series. Generally the best approx- 
imations for the sum result from truncating the series where the terms take the smallest 
absolute values (see e.g. lfT8l ). Since we work in Matlab with a precision of 16 digits and 
with values of \(\ > 10, we typically consider up to 10 terms in the series. In this case the 
terms corresponding to a w are of the order of machine precision (10~ 14 and below). 

Thus we have constructed approximations to the numerical values of the tritronquee 
solution for large values of \(\. These can be used as in [|24ll to set up an initial value problem 
for the P-I equation and to solve this with a standard ODE solver. In fact the approach works 
well on the real axis starting from positive values until one reaches the first singularity on 
the negative real axis. It is straightforward to check the results of [[24]] with e.g. ode45, the 
Runge-Kutta solver in Matlab corresponding to the Maple solver used in ll24l . If one solves 
P-I on a line that avoids the sector | arg (\ > An/ '5, one could integrate until one reaches once 
more large values of \(\ for which the asymptotic conditions are known. This would provide 
a control on the numerical accuracy of this so-called shooting approach. Shooting methods 
are problematic if the second solution to the initial value problem has poles as is the case for 
P-I. In this case the numerical errors in the initial data (here due to the asymptotic conditions) 
and in the time integration will lead to a large contribution of the unwanted solution close 
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to its poles which will make the numerical solution useless. It is obvious that P-I has such 



poles from the numerical results in ||24| and the property (5.5 ). In [24 J the task was to locate 
poles in the tritronquee solution, and in this case the shooting approach seems to be the only 
available. Here we are studying, however, the solution on a line in the complex plane where 
we know the asymptotic conditions for the affine parameter tending to ±00. 

Thus we use as in [fT5l the asymptotic conditions on lines avoiding the sector | arg£| > 
47r/5 to set up a boundary value problem for y = ±y , y > 10. The solution in the 
interval [— y , y ] is numerically obtained with a finite difference code based on a collocation 
method. The code bvp4c distributed with Matlab, see ll39l for details, uses cubic polynomials 
in between the collocation points. The P-I equation is rewritten in the form of a first order 
system. With some initial guess (we use Vl = — a/C/6 as the initial guess), the differential 
equation is solved iteratively by linearization. The collocation points (we use up to 10000) 
are dynamically adjusted during the iteration. The iteration is stopped when the equation is 
satisfied at the collocation points with a prescribed relative accuracy, typically 10~ 10 . The 
solution for a = i and b = is shown in Fig. [TJ The values of f2 in between the collocation 



a 




15 20 



Figure 1 : Real (blue) and imaginary part (red) of the tritronquee solution to the Painleve I 
equation for ( = iy. 

points are obtained by interpolation via the cubic polynomials in terms of which the solution 
has been constructed. This interpolation leads to a loss in accuracy of roughly one order of 
magnitude with respect to the precision at the collocation points. To test this we determine 
the numerical solution via bvp4c for P-I on Chebychev collocation points and check the 
accuracy with which the equation is satisfied via Chebychev differentiation, see e.g. [|45l . 
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It is found that the numerical solution with a relative tolerance of 10 on the collocation 
points satisfies the ODE to roughly the same order except at the boundary points where it 
is of the order 10~ 8 , see Fig. [2] where we show the residual A by plugging the numerical 
solution into the differential equation for the above example. It is straightforward to achieve 
a prescribed accuracy by requiring a certain value for the relative tolerance. Notice that we 
are not interested in a high precision solution of P-I here, but in a comparison of solutions 
to the NLS equation close to the point of gradient catastrophe of the semiclassical system 
with an asymptotic solution in terms of P-I transcendents. For this purpose an accuracy of 
the solution of the order of 10~ 4 will be sufficient in all studied cases. 



-8 - 

-9 - 
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Figure 2: Error in the solution of the Painleve I equation. 

The quality of the used boundary conditions via the asymptotic behavior can be checked 
by computing the solution for different values of y . One finds that the difference between 
the asymptotic square root and the tritronquee solution is only visible near the origin, see 
Fig. [3] For large x it can be seen that the difference between the square root asymptotics 
and the tritronquee solution reaches quickly values below the aimed at threshold of 10~ 4 . 
It is interesting to note that this difference is actually smaller than the difference between 
the tritronquee solution and the truncated formal asymptotic series except at the boundary, 
where the latter condition is implemented (see Fig. [3]). 

The dominant behavior of the square root changes if one approaches the critical lines 
a = exp(47n/5), b = 0. As can be seen from Fig.[4j the solution shows oscillations on top 
of the square root. The closer one comes to the critical lines, the slower is the fall off of the 
amplitude of the oscillations. We conjecture that these oscillations will have on the critical 
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Figure 3: The plot on the left side shows the absolute value of the difference between the 
tritronquee solution and the asymptotic condition — a/C/6 for a = i and 6 = 0. The plot on 
the right side shows in blue the same difference for x > 10 and in red the difference between 
the tritronquee solution and the truncated asymptotic series. 

lines only a slow algebraic fall off towards infinity. 

The above approach thus allows the computation of the tritronquee solution for a line 
avoiding the sector | arg£| > 47r/5 with high accuracy. The picture one obtains by com- 
puting Vt along several such lines is that there are indeed no singularities in the sector 
| arg^| < 47r/5, and that the square root behavior is followed for large To obtain a 
more complete picture, we compute the tritronquee solution for | arg£| < 47r/5 — 0.05 and 
|C| < R ( we choose R = 20). The boundary data for \(\ = R follow as before from the 
truncated asymptotic series, the data for arg ( = ±47r/5 — 0.05 are obtained by computing 
the tritronquee solution on the respective lines as above. 

To solve the resulting boundary value problem for the P-I equation is, however, computa- 
tionally expensive since we have to solve an equation in 2 real dimensions iteratively. Since 
the solution we want to construct is holomorphic there, we can instead solve the harmonicity 
condition (the two dimensional Laplace equation) for the given boundary conditions. To this 
end we introduce polar coordinates r, and use a spectral grid as described in [45 J: the main 
point is a doubling of the interval r G [0, R] to [-R, R] to allow for a better distribution of 
the Chebychev collocation points. Since we work with values of < <p , we cannot use the 
usual Fourier series approach for the azimuthal coordinate. Instead we use again a Cheby- 
chev collocation method. The found solution in the considered domain is shown in Fig. [5] 
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Figure 4: Real (blue) and imaginary part (red) of the tritronquee solution close to the critical 
line (for a = exp(i(4.7r/5 — .05))) with oscillations of slowly decreasing amplitude. 

and Fig. [6j The quality of the solution can be tested by plugging the found solution to the 
Laplace equation into the P-I equation. Due to the low resolution and problems at the bound- 
ary, the accuracy is considerably lower in the two dimensional case than on the lines. This 
is, however, not a problem since we need only the one dimensional solutions for quantitative 
comparisons with NLS solutions. The two dimensional solutions give nonetheless strong 
numerical evidence for the conjecture that the tritronquee solution has globally no poles in 
the sector | arg£| < An/5. 

7 Critical behavior in NLS and the tritronquee solution of 
P-I: numerical results 

In this section we will compare the numerical solution of the focusing NLS equation for 
two examples of initial data for values of e between 0.1 and 0.025 with the asymptotic so- 
lutions discussed in the previous sections, the semiclassical solution up to the breakup and 
the tritronquee solution to the Painleve I equation. The numerical approach to solve the NLS 
equation is discussed in detail in [|29ll . For values of e below 0.04 we have to use Krasny 
filtering 11301 (Fourier coefficients with an absolute value below 10~ 13 are put equal to zero 
to avoid the excitation of unstable modes). With double precision arithmetic we could thus 
reach e = 0.025, but could not go below. 
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Figure 5: Real part of the tritronquee solution in the sector r < 20 and \<j)\ < 4tt/5 — 0.05. 



7.1 Initial data 



We consider initial data where u(x,0) has a single positive hump, and where v(x, 0) is 
monotonously decreasing. For initial data of the form u(x,0) = A 2 (x) and v(x, 0) = 
where the function A(x) is analytic with a single positive hump with maximum value A , 



the semiclassical solution of NLS follows from ( 2. 1 1 ) with f(u, v) given by 

\ 



; v+iWu 



/(«,«) = 29 



^? p(v)\ (v + o 17 ) 2 + u 



iA 



(7.1) 



/ 



where 



p(v) 



7T 



a;+(?7) 



and where x ± (r]) are defined by A(x ± (r])) = ir\. The formula (7.1 ) follows from results by 
S.Kamvissis, K.McLaughlin and P.Miller in [|26l . 

From f(u, v) it is straightforward to recover the initial data from the equations 



x 



fui fv 



0. 



(7.2) 



Numerically we study the critical behavior of two classes of initial data, one symmetric with 
respect to x which were used in ll34l . and initial data without symmetry with respect to x 
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Figure 6: Imaginary part of the tritronquee solution in the sector r < 20 and \(f)\ < An/ 5 — 
0.05. 

which are built from the initial data studied in 11431 . For the former class the corresponding 
exact solution of focusing NLS is known in terms of a determinant. Nonetheless we integrate 
the NLS equation for these initial data numerically since this approach is not limited to 
special cases, but can be used for general smooth Schwartzian initial data as in the latter 
case. 



7.1.1 Symmetric initial data 

We consider the particular class of initial data data 

u(x,t = 0) = Alsecb 2 x, v (x, t = 0) = — ^tanhx, // > 0. 
Introducing the quantity 



we find that the semiclassical solution for these initial data follows from (2.1 1 ) with 



f{u, v) =| V -\{v- 2M)A + - i(v + 2M)A_ - l -u\ogu 



+ 2 M lo S 



\v + M + A+)(--v -M + A_) 



(7.3) 



(7.4) 
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where 



i 

A± = ( {-K)±Mf + u 



For = we recover the Satsuma-Yajima 11371 initial data that were studied numerically in 
11341 . The function f(u, v) takes the form 



f(u,v) = 3l 



-- + iA ) \/u 



«,-.V, , (-f + ^o) 2 + V / (-f + ^o) 2 + 

(7.5) 

(7.6) 



which can also be recovered from ( |7.1| ) by setting p = i. The critical point is given by 

1 



u = 2A 2 , v = 0, x = 0, t 



2Ar 



Furthermore we have 



f = f 

J uuu > J uu 



AAV 



r = 4Al, ip = 0, 



(7.7) 



where r, ^ are defined in ( |2.13| ). For A = 1 the initial data ( |7.3| ) coincides with the one 
studied in Il43l by A.Tovbis, S.Venakides and X.Zhou. In the particular case ji — 2, Aq — 1 



the function /(tt, u) in (7.1 ) simplifies to 



v 1 

f(u, v) = v - - d -v 2 + u + u log 



In this case the critical point is given by 

V = 0, Uq = 2 + fJL, t 

Furthermore, 

f° =0 f° 



(7.8) 



2 + /x 



x = 0. 



r = , ip — 0. 



7.1.2 Non-symmetric initial data 

Recall that we are interested here in Cauchy data in the Schwartz class of rapidly decreasing 
functions. The above initial data are symmetric with respect to x, u is an even and v an odd 
function in x. To obtain a situation which is manifestly not symmetric, we use the fact that if 
/ is a solution to (2.5), the same holds for derivatives and anti-derivatives of / with respect 
to v and for any linear combination of those. If /„ is an even function in v, this will obviously 
not be the case for a linear combination of / and /„. 
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As a specific example, we consider the linear combination 



f = fi + af 2 , a = const, 



where f\ coincides with (jVT8]) and 

[i 2 (i \ 3/2 

f 2 = 2 u J -v 2 + u - - I -v 2 + uj + u v lo£ 



(7.9) 



The function f 2 is obtained from the integration f 2jV — fi — v. The critical point is given in 
this case by 

. 2 , 1, 1 + 4q 1 a _ 1 + 4a 

w = 4(1 -16a), v = -16a, x = - log - _ — , t = - - - log - _ ; 

thus we have |a| < 1/4. Furthermore, 

o 4a 2 - 1/8 a 

J UUV , rz TT, n 1 J u 



Vl-^a 2 ' Juuu 4Vl-16a 2 ' 
such that 

ay/1 — 16a 2 

r = 8u , i> = - arctan 1/g _ 4q2 ■ 

We determine the initial data corresponding to / for a given value of |a| < 1/4 by 
solving ( 1 .2) for u, v in dependence of x. This is done numerically by using the algorithm of 



0T1 which is implemented as the function fminsearch in Matlab. The algorithm provides an 
iterative approach which converges in our case rapidly if the starting values are close enough 
to the solution, which is achieved by choosing u and v corresponding to fo as an initial guess. 
For a close to 1/4 we observe numerically a steepening of the initial pulse which will lead 
to a shock front in the limit a — > 1/4. For the computations presented here, we consider the 
case a = 0.1 which leads to the initial data shown in Fig. [7J 

The initial data are computed in the way described above to the order of the Krasny filter 
on the interval x 6 [—15, 11] on Chebychev collocation points. Standard interpolation via 
Chebychev polynomials is then used to interpolate the resulting data to a Fourier grid. To 
avoid a Gibbs phenomenon at the interval ends due to the non-periodicity of the data, we use 
a Fourier grid on the interval [— 10n, 10n] to ensure that the function u takes values of the 
order of the Krasny filter. For x < — 15 and x > 11, the function u is exponentially small 
which implies the zero-finding algorithm will no longer provide the needed precision. Thus 
we determine the exponential tails of the solution to leading order analytically. We find for 
x — > — oo 



u = v 2 exp 



2(x — av + 
av + + 1 



. 2(x — av+)\ ( , , 2(x — av + ) . 

v + - v + exp ^ -±l 2 log K + ^ — ±1 , (7.10) 

( av + + 1 J \ av + + 1 ' 



26 



1.5 




x 

2 i 1 1 r 




-3 ' 1 1 1 1 1 1 1 

-4-3-2-1 1 2 3 4 



Figure 7: Initial data for the NLS equations without symmetry with respect to x. 



and for x — > +oo 

o / 2(x + a) 

u = v_ exp 

\ av- + 1 

/ 2(x + a)\ ( , N 2(x + a)\ 
v = v--v- exp - 1 — ( 2 log - 1 ; , (7.11) 

\ at)_ + 1 J \ av- + 1 / 

where v± = ± a — l)/a. The initial data for the NLS equation in the form \^ = 
yfu exp(iS/e) are then found by integrating v on the Chebychev grid by standard integration 
of Chebychev polynomials. The exponential tails for S follow from (7.10) and ( |7.1 1[ ). The 



matching of the tails to the Chebychev interpolant is not smooth and leads to a small Gibbs 
phenomenon. The Fourier coefficients decrease, however, to the order of the Krasny filter 
which is sufficient for our purposes. Thus we obtain the non-symmetric initial data with 
roughly the same precision as the analytic symmetric data. 



7.2 Semiclassical solution 

For times t t c , the semiclassical solution gives a very accurate asymptotic description for 
the NLS solution. The situation is similar to the Hopf and the KdV equation JT9). We find 
for the symmetric initial data for t = t c /2 that the L ra norm of the difference between the 
solutions decreases as e 2 . More precisely a linear regression analysis in the case of symmetric 
initial data (for the values e = 0.03, 0.04, . . . , 0.1) for the logarithm of this norm leads to an 
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error proportional to e a with a = 1.94, a correlation coefficient r = 0.9995 and standard error 
a a = 0.03. In the non-symmetric case, we find a = 1.98, r = 0.999996 and a a = 0.003. 

Close to the critical time the semiclassical solution only provides a satisfactory descrip- 
tion of the NLS solution for large values of \x — x c \. In the breakup region it fails to be 
accurate since it develops a cusp at x c whereas the NLS solution stays smooth. This behavior 
can be well seen in Fig. [8] for the symmetric initial data. The largest difference between the 
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Figure 8: The blue line is the function u of the solution to the focusing NLS equation for 
the initial data u(x,0) = 2sechx and e = 0.04 at the critical time, and the red line is 



the corresponding semiclassical solution given by formulas (2.4). The green line gives the 



multiscales solution via the tritronquee solution of the Painleve I equation. 

semiclassical and the NLS solution is always at the critical point. We find that the norm 
of the difference scales roughly as e 2//5 as suggested by the Main Conjecture. More precisely 
we find a scaling proportional to e a with a = 0.38 and r = 0.999997 and a a = 4.2 * 10~ 4 . 
For the non-symmetric initial data, we find a = 0.36, r = 0.9999 and a a = 0.002. The 
corresponding plot for u can be seen in Fig. [9} 

The function v for the same situation as in Fig.[8]is shown in Fig. 10 It can be seen that 



the semiclassical solution is again a satisfactory description for \x — x c \ large, but fails to 
be accurate close to the breakup point. The phase for the non-symmetric initial data can be 



seen in Fig.[TT] In the following we will always study the scaling for the function u without 
further notice. 
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Figure 9: The blue line is the function u of the solution to the focusing NLS equation for 
the non-symmetric initial data and e = 0.04 at the critical time, and the red line is the corre- 



sponding semiclassical solution given by formulas (2.4). The green line gives the multiscales 
solution via the tritronquee solution of the Painleve I equation. 



7.3 Multiscales solution 



It can be seen in Fig. [8] and Fig. 10 that the multiscales solution (5.21 ) in terms of the 
tritronquee solution to the Painleve I equation gives a much better asymptotic description 
to the NLS solution at breakup close to the breakup point than the semiclassical solution 
for the symmetric initial data. For larger values of \x — x c \, the semiclassical solution pro- 



vides, however, the better approximation. The rescaling of the coordinates in (5.21 ) sug- 
gests to consider the difference between the NLS and the multiscales solution in an interval 
[— 7e 4//5 , 7e 4//5 ] (we choose here 7 = 1, but within numerical accuracy the result does not 



depend on varying 7 around this value). These intervals can be seen in Fig. 12 We find that 
the Lqo norm of the difference between these solutions in this interval scales roughly like 
e 4 / 5 . More precisely we have a scaling e a with a = 0.76 (r = 0.998 and a a = 0.019). 

For the non-symmetric initial data, the situation at the critical point can be seen in Fig. [9] 



and Fig. 11 Again the multiscales solution (5.21 ) gives a much better description close to 
the critical point than the semiclassical solution. However, the approximation is here much 
better on the side with weak slope for u than on the side with strong slope. We consider 
again the Loo-norm of the difference between the multiscales and the NLS solution in the 
interval [— 7e 4//5 , 7e 4 / 5 ]. The scaling behavior of the solution can be seen in Fig. 13 For 
7 = 1 we find a = 0.71, r = 0.998 and a a = 0.02. These values do not change much for 
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Figure 10: The blue line is the function v of the solution to the focusing NLS equation 
for the initial data w (x) = 2 sechx and e = 0.04 at the critical time, and the red line is 
the corresponding semiclassical solution given by formulas (2.4). The green line gives the 
multiscales solution via the tritronquee solution of the Painleve I equation. 



larger 7. For smaller 7 there are not enough points to provide a valid statistics. The value 
of a smaller than the predicted 4/5 is seemingly due to the strong asymmetry in the quality 
of the approximation of NLS by the multiscales solution as can be seen from Fig. [9} In the 
considered interval, the deviation is already so big that the scaling no longer holds as in the 
symmetric case. To study the scaling with a reliable statistics would, however, require the 
use of a considerably higher resolution which would be computationally too expensive. 

Going beyond the critical time, one finds that the real part of the NLS solution continues 
to grow before the central hump breaks up into several humps. Notice that the multiscales 
solution always leads to a function u that is smaller than the corresponding function of the 
NLS solution at breakup and before. This changes for times after the breakup as can be 



inferred from Fig. 14 which shows the time dependence of the NLS and the corresponding 
multiscales solution for the non-symmetric initial data. The approximation is always best at 
the critical time. 



To study the quality of the approximation (5.21 ), we use rescaled times. The scaling of 



the coordinates in (5.21 ) suggests to consider the NLS solution close to breakup at the times 
t±(e) with 

*±(e) =t c + u /r - J(uo/r) 2 ±e*/ 5 l3, (7.12) 
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Figure 11: The blue line is the function v of the solution to the focusing NLS equation 
for the non-symmetric initial data and e = 0.04 at the critical time, and the red line is 
the corresponding semiclassical solution given by formulas (2.4). The green line gives the 
multiscales solution via the tritronquee solution of the Painleve I equation. 



where (5 is a constant (we consider j3 = 0.1). We will only study the symmetric initial data in 



this context. Before breakup we obtain the situation shown in Fig. 15 It can be seen that the 
multiscales solution always provides a better description close to x c than the semiclassical 
solution, and that the quality improves in this respect with decreasing e. We find that the 
norm of the difference scales in this case as e a with a = 0.55 (r = 0.994 and a a = 0.03). 



The situation for times after breakup can be inferred from Fig. 16 Close to the central 
region the multiscales solution shows a clear difference to the NLS solution. But it is inter- 
esting to note that the ripples next to the central hump are well approximated by the Painleve 
I solution. The L ro norm of the difference between the two solutions scales roughly like e. 
More precisely we find a scaling e a with a = 1.02 (r = 0.9999 and a a = 7.7 * 10~ 3 ). 



8 Concluding remarks 

In this paper we have started the study of the critical behavior of generic solutions of the fo- 
cusing nonlinear Schrodinger equation. We have formulated the conjectural analytic descrip- 
tion of this behavior in terms of the tritronquee solution to the Painleve-I equation restricted 
to certain lines in the complex plane. We provided analytical as well as numerical evidence 
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Figure 12: The blue line is the solution to the focusing NLS equation for the initial data 
uo(x) = 2 sechx at the critical time, and the green line gives the multiscales solution via the 
tritronquee solution of the Painleve I equation. The plots are shown for two values of e at the 
critical time. 

supporting our conjecture. In subsequent publications we plan to further study the Main 
Conjecture of the present paper by applying techniques based, first of all, on the Riemann - 
Hilbert problem method 11261 l43l l44l and the theory of Whitham equations (see |[T9ll for the 
numerical implementation of the Whitham procedure in the analysis of oscillatory behavior 
of solutions to the KdV equations). The latter will also be applied to the asymptotic descrip- 
tion of solutions inside the oscillatory zone. Furthermore we plan to study the possibility 
of extending the Main Conjecture to the critical behavior of solutions to the Hamiltonian 
perturbations of more general first order quasilinear systems of elliptic type. Last but not 
least, it would be of interest to study the distribution of poles of the tritronquee solution in 
the sector | arg £| > ^ and to compare these poles with the peaks of solutions to NLS inside 
the oscillatory zone. The elliptic asymptotics obtained by Kitaev 11281 might be useful for 
studying these poles for large 

In this paper we did not study the behaviour of solutions to NLS near the boundary u = 0. 
Such a study is postponed for a subsequent publication. 
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Figure 13: The blue line is the solution to the focusing NLS equation for the non-symmetric 
initial data at the critical time, and the green line gives the multiscales solution via the 
tritronquee solution of the Painleve I equation. The plots are shown for two values of e 
at the critical time. 
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